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Abstract. On the basis of the Lorentz equations of motion, the orbit of a charge 
driven by a generic E.M. field with planar symmetry is formulated and analyzed within 
the framework of a Lorentzian geometry with a curvature whose order of magnitude 
i-^ . is parametrized by the radiation intensity and frequency. This reformulation leads to 

(i) a demonstration of the integrability of the particle motion in a plane-wave field 
1^ ' as a result of the vanishing of the curvature, (ii) a manifestation of the parametric 

, dependence of the dynamical response of the particle orbit to the E.M. field on the 

impulse factor in terms of local stability and the occurrence of parametric resonance 
and (iii) a mathematically precise meaning to the ponderomotive oscillation center 
of the charge executing oscillatory motion in a sufficiently low impulsive E.M. field, 
which is subsequently used to raise (iv) a discussion of the domain of applicability of 
the ponderomotive approximation. 

OO 
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. . 1. Introduction 

I The starting point for many investigations on nonlinear radiation-matter interaction has 
been the orbit of a single charged particle in an E.M. plane wave studied on the basis of 
the Lorentz equations of motion or its equivalent Lagrangian formulation. For example, 
see [H 121 [3l m |5l [6] . The primary aims of this paper are (1) to give an exact differential 
geometric formulation of the Lorentz equations of motion in an E.M. field with planar 
symmetry and (2) to analyze the dynamics of the orbits of charged particles in response 
to such an E.M. radiation within the geometrical framework. It is worth noting that the 
class of E.M. fields with planar symmetry includes arbitrary plane-wave fields, standing- 
wave fields, and E.M. fields of two-counter-propagating waves as particular instances. 
See (13.11) . (14. 3p . and (15.11) below. It is shown that the complexity of particle dynamics 
can he encapsulated economically into a curvature function and associated geometrical 
quantities. To demonstrate this conceptual novelty, we show that the integrability of 
the particle motion in an arbitrary E.M. with planar symmetry is a consequence of 
the vanishing of the curvature; the widely known admission of analytic solutions by 
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the Lorentz equations for arbitrary plane-wave fields is therefore a manifestation of the 
flatness of a manifold. 

Another demonstration of the conceptual novelty, which is more important, is 
the availability of the Jacobi equations for geodesic flows and the associated curvture 
features. This way the stability question of particle orbits is translated into the measure 
of the strength of the geodesic spread. To illustrate this, we consider the stability 
of particle orbits in a standing wave field by analyzing the associated Jacobi fields 
and indicate the interrelation of the (in)stability of particle orbits and (non)occurene 
of parametric resonance. Finally, as an important application of the differential 
geometric framework, we employ a mathematical averaging of the Jacobi equations 
and the presence of parametric resonance to discuss the domain of applicability of the 
ponderomotive oscillation center dynamics of charged particles. 

The outline of this paper is as follows. In Sec. Ej the identification of the physical 
orbits of a natural mechanical system with the geodesies in the configuration space 
with a suitable semi-Riemannian metric is established. In Sec. [3l the reduction of the 
interaction of a charged particle with a generic planar symmetic E.M. field to a natural 
mechanical system is made. This sets up the foundation of the differential geometric 
analysis of the radiation-particle interaction in this paper. 

In Sec. m the well-known integrability of the particle motion in a plane-wave field 
is shown as a result of the vanishing of a curvature function. In addition, the analytic 
solutions of the Lorentz equations for an arbitrary plane-wave field is obtained from the 
corresponding geodesic equations. 

In Sec. [5l the stability question of the particle motion in an E.M. field of standing 
wave is studied. The effect of the dynamical response of a particle of charge q and mass 
m placed in an E.M. field of radiation frequency u and field amplitude S is very often 
expressed by the magnitude of the dimensionless impulse factor 



which is usually also called the electron quiver velocity if the charge is an electron or 
the intensity parameter [71 [8]. There is no exception when the particle interacts with 
a generic E.M. field with planar symmetry. However, we will show in this paper that 
the stability of the orbit of a particle interacting with the E.M. field is not completely 
determined by the order of magnitude of the impulse factor t], but by the different 
numerical subintervals within which the value ofrj lies. In particular, a highly impulsive 
E.M. field (77 ^ 1) may not lead to instability of orbits of particles. This, as we will 
show, is a consequence of the (non)occurrence of parametric resonance in the system 
consisting of the particles and the E.M. field, where the impulse factor 77 is the parameter 
of the system that determines the (non)occurence of such a parametric resonance. 

We also employ the Jacobi equations to consider the applicability of the oscillation 
center dynamics of charged particles. The oscillation center dynamics of a charged 
particle in an E.M. field, governed by the so-called ponderomotive force, has been 
the basis of many investigations of the particle motion in response to the field 
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[lElIIllllIinillllinillSlIIlE]. As another significant application of the 
differential geometric framework, we provide a mathematically precise meaning to the 
ponderomotive oscillation center of a charged particle executing oscillatory motion in a 
sufficiently low impulsive E.M. field (?7 <^ 1). This we do by indicating the quantitative 
relation of a Jacobi field and the associated ponderomotive oscillation center to the 
corresponding solution of a mathematically precise averaged version of the Jacobi 
equations. With the aid of this setup, it is shown that when a parametric resonance 
occurs, the ponderomotive approximation breaks down. 

Finally, in Sec. [6l we summarize the results of this paper. 



2. Semi-Riemannian Formulation of Natural Mechanical Systems 

This is a preparatory section in which we derive the identification of the physical orbits of 
a natural mechanical system, a dynamical system that can be described by a Lagrangian 
of the form 

^ V rfr j " 2 dr dr ^ ^ ^ 

where kAB is a metric based kinetic energy tensor and $ is a potential, all of which 
depend only on the coordinates x"^, with the geodesies in the configuration space 
endowed with a suitable semi-Riemannian metric. We will also recall the geodesic 
equations and the Jacobi equations for geodesic spread for later reference. As a 
convention, all capital Latin letters A, B, C, and so on, whenever appearing as subscripts 
or superscripts, always run over the values 0, 1, ■ ■ ■ , — 1, where is the number of 
degrees of freedom of the original dynamical system. 

As is well-known, a natural mechanical system can also be described by the 
Hamiltonian 

nipA,x^) = ^^^ + ^x% (2.2) 

where pa = dC/ d{dx^ / dr) = kAsdx^ /dr are the generalized momenta, and the 
Hamiltonian is an integral of motion. We are particularly interested in the dynamics of 
timelike orbits along which kAB [dx'^/dr) {^dx^ /dr) = if— $ < for a fixed Hamiltonian 
value = H in the case kAB is Lorentzian. In this case, it can be shown that the timelike 
orbits are geodesies in the configuration space endowed with the metric tensor 

gAB = -2{H - ^)kAB, (2.3) 

called the Jacobi metric. This is a classical result for the Riemannian case, a derivation 
based on Hamiltonian's least action principle of which can be found in [151 [HI 113 118] . 
However, we show that it comes as a consequence of the duality of the Hamiltonian T-L 
and the orbit parameter r. Explicitly, we consider r as an additional degree of freedom 
and reparametrizes the orbits by a new parameter A with the property that dr/dX > 0. 
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Since J^I^Xo) ^^'^ = ^ {dr/dX) dX, the Lagrange's equations in the x'^-directions for 
the extended Lagrangian 

2;^ ^ — ^ = £ (x' ^"^^ 



d\ dX J \ dX dX J dX 4^ 

ksc dx^ dx'^ dX ^, A.dr 

= $1 T I — 

2 dX dX dr ^ ' dX 

are equivalent to the ones for the original Lagrangian C Then since r is cyclic with 

respect to C^, 

dC' kAEdx^dx"" fdXV _^ ^2 5) 



d{dT/dX) 2 dX dX \dT J 

is an integral of motion. In fact, comparing (12. 2p and ( 12. 5p . we see that has the value 
—H, so that 

kABdx^dx'' fd_xV ^^^^ (2.6) 



2 dX dX \dT 



Solving for 



dr / kAB dx"^ dx^ 



dX y 2{H - $) dX dX 
from (12. 6 p and applying the classical Routh's procedure yields the Routhian 

CUx^,^]= c 



(2.7) 



dx'^\ ^„ dC dr 



dX J d{dT/dX) dX 



(2. 



dX dX 

that acts also as a Lagrangian for timelike orbits of the original dynamical system. This 
implies that the timelike orbits are geodesies of the configuration space endowed with the 
Jacobi metric l{2.3\) . 

In local coordinates, the geodesies equations are given by 

d'^x^ dx^ dx^ , , 

where 



is the usual affine parameter and 



A 

^ Bc — ~^{9db,c + 9dc,b — Qbcd) 

(2.11) 

are the Christoffel symbols. As is widely recognized, the local stability of a fiducial 
geodesic is governed by the Jacobi separation field J"^ that evolves according to the 
Jacobi equations 

V72 jA J B J D 
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where 



ds 



dJ"^ 
ds 



+ r 



dx^ 



BC 



ds 



-J 



c 



is the covariant derivative of and 



-pA _ -pA I -pA -pE 



ED'- BC 



(2.13) 



(2.14) 



BCD — BD,C BC,D "r EC'' BD 

are the components of the Riemann curvature tensor. More exphcitly, the Jacobi 
equations are given by 

d'^J'^ r.^A dx^ dJ^ ^4 dx^ dx^ 



-J 



D 



0. 



(2.15) 



' " ds ds ^^'^ ds ds 
This way we have translated the stabihty problem of orbits in the original mechanical 
system into geometric language. 

Specializing to the case k^B = Vab, so that qab = ^^"Vab, where 
In [2($ - H)] 



a 



(2.16) 



and 



AB 



-10 
1 
1 














1 
1 



(2.17) 



is the usual Lorentzian metric tensor, we have 



r^BC,D — ^^b'^,CD + ^^c^,BD — VBC<^'^.D 



and 



(2.18) 
(2.19) 

(2.20) 

ds"^ ' ds ds '^^ ds ds ^ ^ 

and, upon using the relation ds/dr = e^" that follows from ^M), fl2A0|) . and (l216|) . by 



R\cD = ^%{(^,BC - <^,B<y,c) - 5%{<^,BD - Cr,BCr,D) + (S^dVBC " '^^c'7Bo)^''',i? 

So the geodesic equations are given by 



d'^x^ da dx^ a dx^ dx'^ 

+ 2- ; a'^^VBC— — = 



cr'^TlBc- 



while parametrized by the original orbit parameter r. Using ( 12. 6 p with A = r and 
kAB = Vab, it can be shown that fl2.22p is equivalent to 

d'^x^ 



+ a'^e^'^ = 0. 



(2.23) 
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Since 

cT,Ae'"' = (2.24) 
by fl2.16p . f l2.23p is in turn equivalent to 

r]^^^^B = 0. (2.25) 



These are simply the Euler's equations associated with the original Lagrangian C with 
kAB = Vab- 

On the other hand, since e^^rjAB {dx^/ds) [dx^/ds] = -1 by (Q, (IZIO]), and 
f l2.16p . the Jacobi equations for timehke orbits are given by 

d'^J^ dadJ^ dx^ d , ^p,. . dx^ dJ'^ a / x 

^ + ^TsltT + ^ItTTs ("^--^ ) - '"^^-^^^ + a'^J-e-^ = 0. (2.26) 

Using the relation ds/dr = e^'^, the Jacobi equations are also given by 

d^J^ fdx"^ A dx^\ dJ^ ( dx^ dx^ A 2a\ .R , 

Denoting the coordinate x^ by t, we call t the lab time and an orbit along which x^ 
remains constant for all A 7^ a spatially constant orbit. On a spatially constant orbit, 
a^A vanishes identically for all A 7^ by any of the geodesic equations f l2.2ip . f l2.22p . or 
f l2.23p . and dt/dr = e" by (12. 6p and (I2.16p . Hence on such an orbit, the Jacobi equations 
can be reduced to 

d^J^ dadJ^ 



+ g-2<x _ (2.28) 

ds ds 

+ e'^a = 0, (2.29) 



ds"^ ds ds 
d^J^ 



and 



dr^ 

d^J^ dadJ^ 



while being parametrized by the affine parameter s, the original orbit parameter r, and 
the lab time t respectively. 



3. Dynamics of a Charge in a Generic E.M. Wave with Planar Symmetry 

We begin this section with the reduction of the interaction of a charged particle and a 
generic E.M. field with planar symmetry to its longitudinal degrees of freedom. Then we 
show that the motion of the particle in the direction longitudinal to the field propagation 
is a natural mechanical system and hence the corresponding orbits are the geodesies in 
the reduced configuration space with a suitable Jacobi metric. 

A generic E.M. with planar symmetry is defined by the vanishing of the longitudinal 
components of its vector potential and its sole dependence on the longitudinal 
coordinates: 



— -^x — 0, Ay — y4j,(t,x), Az — ^2(^,0;). 



(3.1) 
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Together with the Lorentz equations of motion 
d'^x'' „,« fOA^ dA,3\ dx"/ 



= qv'"' ^ - ^ Py-, (3.2) 

dr^ V dxl^ dx'f J dr ^ ^ 

the vector potential determines completely the dynamics of a particle of mass m 
and charge q in the E.M. field. As a convention, all Greek letters a,/3,7, and so on, 
whenever appearing as subscripts or superscripts, always run over the values 0, 1, 2, and 
3. We will also use capital and lower-case Latin letters to denote the longitudinal and 
transverse degrees of freedom respectively. 

The solutions of the Lorentz equations are the world lines 

(x°(r)) = (x^(r); x»(r)) = (t(r), x(r); y(r), z{t)), (3.3) 

each of which is parametrized by the particle's proper time r. The coordinate t = x^ is 
the lab time in units of light traveling distance and is related to the conventional time 

''conv 

by 

t ctconv- (3*4) 

As is generally known, the Lorentz equations are the Euler-Lagrange's equations of the 
Lagrangian 

with the equivalent Hamiltonian 

= n'-'(P^-lA,)(p,-lA-,) ^ ^^ ^^ 

where pa = dC^'^'' /d{dx"'/dT) = rjajBdx^ /dr + qAa/m are the generalized momenta. As 
a consequence of the synchronization of the laboratory clock and a clock comoving with 
the particle whenever the particle is at rest [dx/dr = dy /dr = dz/dr = 0) in the lab 
frame, it follows that Ti^^^ = — ^. 

The planar symmetry of the vector potential induces the cyclic nature of the 
transverse coordinates y and z with respect to the Lagrangian C^'^\ This leads to 
the fact that the transverse momenta 

Pa=Vab^ + -Aa{t,x) (3.7) 

dr m 

are integrals of motion besides the Hamiltonian "H^^^ Thanks to these integrals of 
motion, we can reduce the number of degrees of freedom of the particle dynamics using 
the classical Routh's procedure. To this end, we solve for 

^=(p,-«^),- (3.8) 

dr \ m / 

from f l3.7p and substitute them into 



^(2) ( ^(4) dC('^ dx^ 



dr J d{dx^/dT) dr g-j 

7]Bc dx^ dx^ 
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where the scalar potential is given by 

2 

to get a Lagrangian C^'^^ for the motion of the particle in the direction longitudinal to 
the field propagation. 

The Routh's reduction indicates that the dynamics of the E.M. field- accelerated 
particle is controlled entirely by the equation 

% - (3.11) 

which are the Euler-Lagrange's equations for the reduced Lagrangian C^'^K By contrast, 
subsequent to the launching of the particle, the transverse degrees of freedom have no 
active effect on the particle dynamics in the tx-plane and are governed by the longitudinal 
dynamics through ([5'. 7| j; that is, the transverse degrees of freedom enter the longitudinal 
dynamics of the particle only as parameters by means of the initial data. 

On the other hand, upon comparing fl3.9p to fl2.ll) . we note that more significant 
is the fact that the longitudinal dynamics is a natural mechanical system with the 
Hamiltonian 

^(^nPA,x^) = ^^^|^ + $(t,x), (3.12) 

where pa = dC^"^^ / d{dx^ / dr) = rjAndx^ / dr are the generalized momenta. Also, along 
any physical orbit, we have 

by (13. 6 p and f l3.10p . In consequence, the longitudinal dynamics can be and will be 
identified with the geodesies in the tx-plane endowed with the metric tensor 

^^^9AB = (1 + 2$) VAB = e^^'vAB, (3.14) 

as is given in the more general case by (12. 3p . 

Being defined on a two-dimensional manifold, the geodesic equations have the 
simple form 

d'^u f du\^ , ^ 

^ + 2.„(-j =0. (3.15a) 

"^" aaJ^V-O (3.156) 



ds"^ ' \ds J 
while parametrized by the affine parameter s, where 

u = t — X and v = t -\- x, (3.16) 

are the null coordinates. The components of the Riemann curvature tensor can be 
reduced to 

^'^R^CD = e'^K {6^aVBD - S^^vbc) , (3.17) 
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where 

K = -e-^^'a'^ (3.18) 
is the Gaussian curvature. Thus the Jacobi equations fl2.12p take the form 

V72 jA J B J D 

^,.^^KiSUso-S^,,,c)'-l^jo'-I^^O. (3.19, 

This way one sees that the stability question of the orbit of a charged particle in response 
to an E.M. field with planar symmetry is encapsulated into the product of two geometrical 
quantities, i.e., the conformal factor e^'^ and the Gaussian curvature K. 

4. Dynamics of a Charge in a Curvature-Free Wave Field with Planar 
Symmetry 

It is clear from (13.181) that the tx-plane endowed with the metric ^"^^Qab is flat if and only 
if (j{u^v) = U{u) + V{v), where U and V are smooth functions of one variable. Under 
this circumstance, initially parallel geodesies preserve their separations [IHl §11.5]. In 
addition, with the vanishing of the Gaussian curvature K., the geodesic equations (I3.15al) 
and ( 13.1561) decouple from each other and are equivalent to 

^fe2^M^^=^fe^v^(.)^Uo. (4.1) 



ds \ ds J ds \ ds J 
Consequently, one has 

/ .«^co„st./.. and const./.. (4.2) 

In summary, this shows that the vanishing of the Gaussian curvature K implies the 
integrability by quadratures of the particle motion. 

A curvature-free instance of physical interest is readily furnished by an elliptically 
polarized plane-wave field propagating along the x-axis. The components of the vector 
potential can be represented by 



S3 S \f\ ^2 
At = Ax = 0, Ay{u) = — cosuu, Az{u) = sinuu, (4.3) 

CO OJ 

where S and u are the field amplitude and frequency respectively, while 5 is a 
polarization parameter such that 5 = 0, ±1 for a linearly polarized wave and 6 = ±l/\/2 
for a circular wave. It can be shown that the equations of motion for a particle in such 
a field can be solved exactly [H [3 IH Ej. Since the vector potential depends only on 
n = t — X, it is readily seen from fl2.16p . f l3.10p . and f l3.18p that the induced tx-plane is 
fiat. Thus the widely known integrability of the particle motion is indeed a consequence 
of the vanishing of the Gaussian curvature K. Furthermore, the tranverse motion can 
be obtained by integrating (13. 7p . while the longitudinal dynamics of the particle can be 
determined with the aid of (14. 2 p by integrating 

J 1 + {Py + ri6 cos uu)"^ + (^Pz + rjV 1 — 5^ sin uu^ du = const. J dv, (4.4) 
where we recall rj = qS/mu is the dimensionless impulse factor of the interaction. 
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5. Dynamics of a Charge in a Standing Wave Field with Planar Symmetry 

In general, the Gaussian curvature K as determined from f l3.18p is nonvanishing. Of 
particular interest is the motion of a charged particle driven by the E.M. field of a 
standing wave |2l [lOl [13] with planar symmetry and polarized linearly along the 2;-axis. 
The components of the vector potential are 

g 

At = = Ay = 0, Az{t,x) = — sinwt sincjx, (5.1) 

to 

where S and u are the field amplitude and frequency respectively. In this case, the 
scalar potential for the longitudinal dynamics is given by 

$(t,x) = ^ + (P, - 7]smutsmux)^] . (5.2) 
2 

Here we recall that 

V = ^ (5.3) 

muj 

is the dimensionless impulse factor of the interaction. The different numerical 
subintervals within which the value of rj lies, as we will see, play a significant role 
in the local stability of the physical orbits of the particle. 

By (13.111) . the longitudinal dynamics is governed by the equations 

— — = (P, — nsmutsmux) cosutsmux, (5.4a) 

dr^ Lo^ ' ' ^ ' 

d'^x rj 

-—r = — (Pz — v sin cot sin cox) sin cot cos ujx. (5.46) 
dr^ CO 

If, in addition, we stipulate that Pz = 0, then there are two classes of spatially constant 
orbits. One class of these spatially constant orbits are those along which sin uox = and 
whose x-coordinates are 

TITT 

x(r) = ±— , n = 0,1,2,- ■■, (5.5) 

CO 

and another class of spatially constant orbits are those along which coswa; = and 
whose x-coordinates are 

a:(r) = ± (^n + ^, n = 0, 1, 2, ■ ■ ■ . (5.6) 

We will use the first class of spatially constant orbits to discuss their local stability 
characterized in terms of the numerical subintervals of the impulse factor 77. 

5.1. Dynamics about the Orbits x = dznn/co 

By setting cr = \ In (2$ + 1) in accordance with (I2.16p and (I3.13p . one finds 



rf 00 cos cot sin cot sin? cox 
1 + Py + rj'^ sm cot sm cox 
Tj"^ CO sin^ cot cos cox sin cox 

^'^ = 1 I P2 I 2-2 . • 2 ' (5-7^) 

i- + Py + f] sm cot sm cox 
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, (1 + P^) cos^ ut sin^ ux — rf sin^ ut sin^ ux 

o.tt = iv^y n ^ P2 I 2-2 . ■ 2 ' (S-'^c) 

(1 + -T^ + r/"^ sm cur sm ux)'' 

^(1 + P^) cos cut sin ut cos sin ux 

(1 + -Tjj + r/^ sm cur sm wxj^ 

„ (1 + P^) sin^ cos^ ux — rf sin^ sin^ ux 

(J,xx = M n ^ P2 I 2-2 . ■ 2 r2 • (5-7e) 

(1 + -Tjj + r/"^ sm sm wxj"^ 

In particular, along the spatially constant orbits x = ±mT/u, = 0, 1, 2, ■ ■ -, one has 

= = and a^tt = (^,tx = 0, 

{r]usmutf (5.8) 



(^.xxif^ 3^) 



It follows from these equations and fl2.30p that the Jacobi equations associated with the 
spatially constant orbits x = ±mT/u, n = 0,1,2, ■■ ■, are 



d^J^ (rjusvuuty 
~W ^ 1 + P2 



= 0, (5.9a) 
r = (5.9&) 



while being parametrized by the lab time t. Due to their simplicity in form, they readily 
admit further analysis to address the dynamical properties of neighboring orbits of the 
spatially constant orbits x = ±n7i/u, n = 0, 1,2, By a neighboring orbit of a 
spatially constant orbit we mean one which begins in a sufficiently small neighborhood 
centered around the initial datum for the particular spatially constant orbit in the phase 
space of initial data. 

Without restricting generality, let Py = and let the fiducial orbit be the spatially 
constant orbit x = 0. On rescaling the lab time using the optical cycle (i.e., 27i/u units 
of lab time), fl5.9&l) becomes 

+ (27rr/ sin 27rT) V^' = 0, (5.10) 



where 



dT^ 

„ UJt 



- 2. ("1) 
is the dimensionless time in units of optical cycles. Since {2'jTriy^ s\.v?2'KTdT = oo, 
(]5.96p is oscillatory in the sense that all its solutions have arbitrarily large zeros [20| 
Theorem 2.4.1]. In other words, the fiducial orbit x = has infinitely many conjugate 
points. Hence we expect that neighboring orbits intersect the fiducial orbit infinitely 
often, in particular, whenever the corresponding solutions of f l5.10p are bounded. The 
distribution of these conjugate points may therefore serve as a physically identifiable 
diagnostic property of the dynamics of the neighboring orbits of the orbit x = as 
estimated by f l5.10p as compared to f l5.4ap and ( 15.4&P 

More pertinent to the stability question of the spatially constant orbit x = is the 
fact that (]5.10p is a Hill equation, so that the parametric dependence on the impulse 
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factor rj of the boundedness of its solutions can be deduced from standard Floquet theory 
as follows. Let and Jf be the solutions of flS.lOp such that 



JX 



1 and 



dJ^ 



J^' = at T = 0. 



(5.12) 



dT dT 

Following [21], the characteristic function (f^t]) of flS.lOp is defined as half the trace of 
its fundamental matrix at a period, i.e.. 



tr 



JX 



JX 
Or. 



dJf/dT dJ^/dT 



(5.13) 



Then it is known that all solutions of fIS.lOp are bounded if and only if |v5(?7)| < 1 
Sec. 2.4]. Since the characteristic function </3(?7) depends generally on the value of 77, 
so does the boundedness of solutions of fl5.10p and hence the stability of the spatially 
constant orbit x = 0. Viewing the impulse factor 77 as a parameter of the system, we 
call any value of 77 for which x = becomes unstable a resonant value and say that 
there is a parametric resonance under such a circumstance. 

In fact, |v5(?7)| — 1 has arbitrarily large zeros on the interval (0, 00). Thus there exist 
infinitely many disjoint open suhintervals (t/;^, 77^4,1) C (0,oo), k = 0, 1,2, ordered 
according to their left endpoints, in which \y^{ri)\ < 1 I21\ Ch. VII, ^1-4, %l-5] and 
for values of rj lying within these intervals, the spatially constant orbit x = is stable. 
These subintervals are called the stability zones; the complementary closed subintervals, 
within the interior of each of which |</3(r/)| > 1, are called the instability zones. As an 
additional remark, it follows from a classical result [23] that if 

"1/2 9 



1 

4 



-r = v 



sm'27iTdT < — , 



(5.14) 



then all solutions of f l5.10p are bounded. Hence the first stability zone must be of the form 



(0, ) for some rj-^^ > 



2V2 



0.900316. Numerically, 7/^ ~ 1.147179. In consequence. 



the orbit x = is stable in a sufficiently low impulsive E.M. field (rj <^1). 

5.2. A Mathematically Precise Averaging Principle for the Jacobi Equation 

In a sufficiently low impulsive E.M. field (7/ ^ 1), the dynamics of the neighboring orbits 
of the spatially constant orbit x = can be analyzed using the method of averaging. 
To this end, we note that fl5.10p is equivalent to the linear periodic system 
dJ 



dT 



2Trr]A{T)J, 



where 



" Jl ' 




J2 





J" 



{2Tiriy^dJ''/dT 



(5.15) 



and the coefficient matrix 



A{T) 



1 
-sin2 27rT 
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13 



It = 



where 



(5.16) 



A 







is the average of A over one of its own cycles or over an optical cycle. The solution of 

^10 



the averaged system, i.e., fl5.16p . with the initial condition J(0) = is readily 

obtained as 

J(T) = e'"''^^ J(0) 

cosv^vrr/T \/2 sin \/27rr/T 1 [ Jiq 1 (5-17) 
— sin -\/27r?7T cosy^vrTyT J20 

showing that the averaged system exhibits simple harmonic motions with the frequency 
r\l^2 oscillations per optical cycle. To consider how good the averaged system (15.161) 
approximates the original system (I5.15p . let 



^1 
E2 



J2 ~ J2 



E20 



J20 — J20 



-E=J~J= 

be the error in estimating J using J with the initial error 

= E{0) = J(0) - J(0) 

Calculating E{T) from ( ]5J[5|l . (^M), and (KWf yields 

E{T) - E{0) = 27rv f [A{S){E{S) - E^ 
Jo 

+{A{S) -A)J{S) + A{S)E{0)] dS. 

In consequence, one finds 

\\E{T) - Em 

<2nr^ [ [\\A{S){E{S)-Em\\ + \\{A{S)-A)JiS)\\ + \\AiS)Em] dS 
Jo 



(5.18) 



(5.19) 



< 



27ir]£ ^J{E^{S)-Ewysm^27^S+iE2{S)-E2oydS 
+ nrj [ cos 4:7tS (Jiq cos \/2Txr]S + J2oV^sin \/2tx7]S 

dS 



(5.20) 



^2A/E2osin^27r5 + E|o 



<27rr/^ \\{E{S)-E{Q)\\ c/^ + Trr/ (^| Jio| + | J20I + 2 ||^;(0)||) T. 
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Then an application of Gronwall's lemma Lemma 1.3.3] yields 

\\E{T) - Em < \M + ^\M + 2\\Em ^^,.,t _ 21) 

This equation implies that if J and J agree initially, then 

J(T) = J(T) + 0(e) (5.22) 
uniformly for < T < ln(l + e)/2Tiri in the sense that 

\\J{T)-J{T)\\ < \M±:/l\Me (5.23) 
for any e > 0. Specifically, 

r{T) = J^{T) + 0{e), (5.24) 

where 

1 d.r 



J'^iT) = J^{0)cosV2nr]T 



sin V2tt7]T (5.25) 

T=0 



V2'iT7] dT 

is the average value of = Jio as obtained from the averaged system (15.161) . i.e., 
= Jio, and the initial data for and have been identified. Hence in a sufficiently 
low impulsive E.M. field (rj <^ 1), a charged particle that is launched from a neighboring 
orbit of the spatially constant orbit x = exhibits an almost simple harmonic motion 
and oscillates about x = with a frequency close to rj/ \/2 oscillations per optical cycle. 

5.3. The Landau Decomposition of the Jacobi Equation 
We can reformulate the Jacobi equation, i.e., f l5.10p . as 

^ = 2'K^rf{-l + cos47rT) (5.26) 

Thus the dynamics of the Jacobi field J^' can be viewed as the motion of a particle 
subject to a time- independent potential 

U{r) = {Tirirf (5.27) 

and a force 

f{r,T) = 2TT'^rf cos AtiT ■ r (5.28) 

which varies in time with the frequency 2 oscillations per optical cycle. The frequency of 
/ is considered high in the sense that it is of higher order of magnitude of the frequency 
of the motion of the particle in the absence of the force field /, i.e., 2 ^ ?7/a/2 in this 
case. Note that this is equivalent to requiring that the E.M. field is of low impulse, i.e., 
7] < 2^2. It is a common practice in plasma physics to decompose the original motion. 



Differential geometric analysis of radiation-particle interaction 



15 



i.e., the Jacobi field in our case, into a secular component X describing the oscillation 
center orbit and an oscillatory component ^ that is nearly periodic [71 [261 1121 E]: 

r{T)=X{T) + aT) 



with 



~dT^ 
dT^ 



-27rV(^ + 2sin2 SttT-^, 
2TT^r]^cos4TTT ■ X. 



(5.29) 



We call the decomposition {X,C,) a Landau decomposition. This is a non-unique 
representation of which depends on a choice of X and ^ satisfying the last two 
coupled equtions of (15.291) . 

Of particular interest is the following exact representation. It is obtained by 
applying the initial conditions J^'(O) and dJ^ /dT\j,^Q to X and dX/dT and having 
the oscillatory component ^ and its velocity d^/dT vanish initially. This gives rise to 
the compatibility condition 

dJ"" 



and 



X{0) 



m 



0, 



d^ 

dT 



T=0 



0. 



dT 



T=0 



(5.30) 



T=0 



They determine an exact solution to the last two coupled equations of (I5.29P and thus 
constitute the exact representation of the Landau decomposition (X, C.) expressed by 
(I5.29p . With the compatibility condition, each Jacobi field .P is uniquely determined by 
a Landau decomposition and vice verse. 



5.4- The Averaging of the Landau Decomposition 

On the assumption that the E.M. field is of low impulse, the Laudau decomposition 
(I5.29P can be analyzed using the method of averaging as in the case of (I5.10p . Together 
with the compatibility condition (15. 30 p . we will show that the Jacobi field and 
the oscillation center X agree on average in a mathematically precise sense, thereby 
quantifying the smallness of the amplitude of oscillatory component ^ as compared 
to that of X. Moreover, the result of the calculation is a preliminary in yielding a 
mathematically solid meaning to the ponderomotive approximation of the oscillation 
center X in a sufficiently low impulsive E.M. field (77 <^ 1). This will be discussed in 
more detail in the next section. 

To initiate the averaging process on the Landau decomposition (I5.29p . we note the 
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equivalence of the last two equations of fl5.29p to the linear periodic system 
dK 



dT 



V2nriB{T)K, 



where 



K 



X 

-1 



[V^irr]) dX/dT 
{V27rriy^ d^/dT 



(5.31) 



and the coefficient matrix 



5(T) 



1 

-1 -2sin^27rr 

1 

cos47rT 

has period half an optical cycle. The averaged system corresponding to ( I5.3ip is 
dK 



dT 



where 



B 



(5.32) 



10 
-10-10 
1 


is the average of B over one optical cycle. The solution of the averaged system with the 

Kio 



initial condition K{0) = — is given by 



K20 



(5.33) 



K{T) 

COS \^nr]T sin y/2nriT cos \/27it]T — 1 sin ^J^ixriT — \f2'nriT 

— sin ^J^-nrfT cos y^TrrjT — sin ^/2nrjT cos y^nrfT — 1 

1 V^TrrjT 

1 

Note that K3q{T) = K^o + \/2TTr]KioT. In order to be consistent with the stipulation 
that the average of ^ vanishes over each optical cycle, we necessitate 

Kso = Km = 0. (5.34) 



X 
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It follows from this that 

Ks{T) = K,{T) = (5.35) 

identically. This renders K purely oscillatory and thus gives a heuristic justification for 
the compatibility condition f l5.30p . With this condition on any solution of the averaged 
system (15.32 p . we analyze how good the averaged system approximates the original 
system (I5.3ip by considering the error 



F, 

with the initial error 



K -K 



K2-K2 



(5.36) 



-^10 

-^20 
F40 



F(0) = K(0) -K(0) 



K20 — K20 





where we have set -^'3(0) = -^'4(0) = in view of the compatibility condition (I5.30p . As 
for (I5.19p . one finds 



F{T) - F(0) = V27r7] [ [B{S){F{S) - F{0)) 

Jo 

+{BiS) - B)K{S) + A{S)F{0)] dS. 



(5.37) 



In view of this, one has 
||F(T)-F(0)|| 

BiS)iFiS) 



Jo 

< V2nr] / [( 
Jo 



{B{S) - B) K{S)\\ + ||fi(5)F(0)||] dS 



F,{S) - Fio)'(l + cos^ AnS) + {F2{S) - F, 



20 J 



+^3(5)2 sin^ 2ttS + Fi{Sf + 2(Fi(5) - ^10)^3(5) sin^ 2tiS] ^'^ dS 

+ V2nr] J I cos AnS (^Kio cos \/2TxrjS + K20 sin \/2'KrjS^ 

+ [Fl{l + cos^ AtiS) + F|o] dS 

<2V2ivri [ \\F{S)-F{0)\\ dS + V27rri{\Kio\ + \K2o\ + \\F{0)\\)T. 
Jo 

Hence an application of Gronwall's lemma [25| Lemma 1.3.3] yields 



(5.38) 



|F(T)-F(0)||< 



20 



|i^(0)|| 



,2V2nriT 



This equation indicates that if K and K agree initially, then 
K{T) = K{T) + 0{e) 



(5.39) 
(5.40) 
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uniformly for < T < ln(l + e)/2\^'7iri in the sense that 



\K- 



10 1 



20 1 



\\K{T) - K{T)\\ < 
for any e > 0. In particular, one concludes that 
X{T) = X{T) + 0{e). 

Here 

1 J V 

X{T) = X{0) cos V27Tr]T 



sm 



(5.41) 
(5.42) 
(5.43) 



T=0 



V27rr/ dT 

is the average value oi X = K\ as obtained from the averaged system fl5.32p . i.e., 
X = Ki, and the initial data for the Landau decomposition (X,^) and its averaged 
version (X,^) have been identified. Thus, combined with the compatibility condition 
fl5.30p . f l5.25p and fl5.43p yield the salient consequence that the original Jacobi field .P 
and the oscillation center X in the Landau decomposition (X, ^) agree in the sense of 
average: 

l^iT) =X{T). 

Hence 

r{T) = XiT) + 0(e), 
m =0{e) 

uniformly for < T < ln(l + e)/2\/2Tiri in the sense that 

1 



(5.44) 

(5.45a) 
(5.456) 



\r{T)-XiT)\ = UT)\< \\.T 



A/27r?7 



dT 



T=0 



(5.46) 



for any e > 0. This quantifies the smallness of the rapid oscillation ^ as compared to 
the oscillation center X. 



5. 5. The Ponderomotive Approximation of the Landau Decomposition 

Based on the compatibility condition (15. 30 p . we will quantify the closeness of the 
oscillation center X in the Landau decomposition (X,^) and its ponderomotive 
approximation Xp, thereby providing a mathematically solid formulation for the 
ponderomotive oscillation center Xp. To this end, we first calculate the ponderomotive 
oscillation center Xp following (TJ [2Sl [121 [S] • More sophisticated derivations that lead 
to equivalent results can be found in [21 [6] . In addition, the result of this calculation 
leads to an understanding of the breakdown of the ponderomotive approximation when 
a parametric resonance occurs. This will be expounded in the next section. 

While regarding X as constant over a given optical cycle and taking into account 
of the zero average of ^ over each optical cycle, we integrate the last equation of (15.291) 
to obtain 

UT) = -^cos47rT-X(T). (5.47) 
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Here we have used the subscript 'p' to distinguish C,p obtained this way from the exact 
rapid small oscillation C, in the Landau decomposition [X,^). Then the evolution 
equation of Xp is obtained by first replacing X and ^ in the second equation of f l5.29p 
by Xp and C,p respectively and then by averaging the resulting equation. The result is 

d^X„ 



2n'^7]'^Xp - Att'Y [ -— cosAttT ) sin^ 27rT ■ Xp 



.2^2 



V 



2 2 

r f] 
7r^?7^(16 + 



-2'iT^ri'Xp -Xp 



-X„ 



(5.48) 



Consequently, 



Xp{T) = X{0) cos \ 2 + ^7r7]T 



1 dX 



sm 



T=0 



-TTrjT. 



(5.49) 



16 + 77^ TTT] dT 

Here the initial data for X and Xp have been identified. Then, in view of ( I5.43P and 
dEH]), one has 



\X{T)-Xp{T)\<2' ■ 
by the mean value theorem. This implies that 
\X{T) - Xp{T)\ < 

and thus 

X{T) = Xp{T) + 0{e 



xm 



1 


dX 




nrj 


dT 


T=oJ 



TiriT 



X 



1 


dX 




TlTj 


dT 


T=oJ 



(5.50) 



(5.51) 



(5.52) 



uniformly for < T < s/27iri for any e > 0. Finally, combining fl5.4ip and fl5.5ip yields 

1 



\X(T)-Xp{T)\< 
showing that 



^ + \/2 + i^ 11^(0)1 + 



2^2 



1 


dX 




TTrj 


dT 


T=0 



(5.53) 



X(T) = Xp{T) + 0(e) (5.54) 

uniformly for < T < min |ln(l + £:)/2A/27r?7, £:/27r?7} ^^Y £ > 0- This gives a 
mathematically rigorous statement about the closeness of the oscillation center X in the 
Landau decomposition {X, C,) to its ponderomotive approximation Xp, whose evolution 
is governed by f l5.48p instead, in a sufficiently low impulsive E.M. field {rj <^ 1). 



5. 6. Breakdown of Averaged Landau Decomposition and Ponderomotive Approximation 
due to Parametric Resonance 

We say that the mathematical averaging principle and the ponderomotive approximation 
for the oscillation center X in the Landau decomposition (X, ^) break down when 
f l5.42p and fl5.54p respectively become inapplicable as T — oo. We will exhibit such 
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a breakdown when a parametric resonance occurs in the original Jacobi field J^. This 
we do by first relating the exponential order of to the ones of X and ^. More explicitly, 
we show that if any of the members of the Landau decomposition {X, ^) has a positive 
exponential ordeii||, then so has the other element the same exponential order. Thus it 
makes sense to talk about the exponential order of the Landau decomposition (X,^) 
instead of the exponential order of any of its elements. Furthermore, the associated 
Jacobi field possesses the same exponential order as the Landau decomposition 
(X,^). On the other hand, when a parametric resonance occurs, the Jacobi field has 
exponential order /i for almost all initial data, where fi is the positive Floquet exponent 
associated with the Jacobi equation, i.e., f lS.lOp . Together with the relation among 
the exponential orders of the Jacobi field and its associated Landau decomposition 
(X,^), this implies that both the oscillation center X and the oscillatory component ^, 
assumed small in the ponderomotive approximation, diverge exponentially in amplitude 
as T ^ CO with exponent no smaller than the Floquet exponent /x for almost all initial 
data in the case when a parametric resonance in takes place. In fact, both dX/dT 
and d^/dT also diverge exponentially in amplitude as T — )■ oo with no smaller exponent 
under the same circumstances. These divergences suggest an inconsistency with f l5.42p 
and f l5.54p as T — )■ oo, where the oscillation center X in the Landau decomposition (X, ^) 
is bounded in both its averaged and ponderomotive approximations. In other words, 
these divergences suggest the breakdown of the mathematical averaging principle and 
the ponderomotive approximation for X. The primary aim of this section is to provide 
a rigorous mathematical justification of these suspicions. 

To begin with, let us relate the exponential orders of the Jacobi field and the 
members in the associated Laudau decomposition (X, 0- To this end, note that the 
evolution equations, i.e., the last two equations of fl5.29p . of the Landau decomposition 
(X, ^) satisfy the following equivalent integral equations 



First of all, note that these equations imply that if X{T) has exponential order a > 
for T > 0, then so have ^{T) and J^{T). To see this, let Mi > be such that 

I If for some a, M e M, |/(T)| < Me"^ for all T > 0, then the function / is said to have the exponential 



X{T) =X{0) cosV27rr/T + 



1 dX 
V27Tr] 'dT 



T=0 





(5.55) 



order a. 
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\X{T)\ < Mie"^ for all T > 0. Then 



\aT)\<m\ + 


d^ 

dT 


TH 

T=0 




= \m\ + 


di 


TH 

T=0 




dT 





-aT-l) 



<\m\ + 



27r'77'Mi 



d^ 



dT 



27rV^i 



T=0 



(5.56) 



27rVMig.T 



< 



ie(o)| + 



di 



dT 



T=0 



by the second equation of fl5.55p . from which it follows that 



\r{T)\ < 



Mi + |e(0)| + 



a 



di 



dT 



+ 



T=0 



a" 



(5.57) 



for all T > 0. 

Secondly, if if C,{T) has exponential order /3 > for T > 0, then so have X{T) and 
J'^iT). Explicitly, let |^(T)| < Mae^^ for T > 0, where Ms > 0. Then 



|X(T)| < \X 



\X 



+ 



< \x 



< 



\x 



1 

1 

1 



dX 



dT 
dX 



2V2' 



dT 
dX 



T=0 



T=0 



dT 
dX 



T=0 



dT 



T=0 



nrj 
2^/2-117] 

2^/21X1] 
4\/27r?7 



dS 



(5.5J 



+ 



and 



\riT) \ < 



M2+ |X(0)| 



y/2TTr] 



dX 



dT 



+ 



T=0 



A^/27^r| 



(5.59) 



for all T > 0. 

We are now ready to employ ( I5.57P and f l5.59p to justify the statement that 
when a parametric resonance occurs in the original Jacobi field , i.e., when the 
characteristic function v'(?7) of f l5.10p has an absolute value |v5(?7)| > 1, both X and 
C, diverge exponentially in amplitude as T — )■ 00 with exponent no smaller than fi for 
almost all initial data for , where fi is the positive Eloquent exponent associated 
with the Jacobi equation Ii5.10\) . As we have pointed out earlier, this suggests the 
breakdown of the mathematical averaging principle and the ponderomotive approximation 
for the oscillation center X because there both the averaged and the ponderomotive 
approximations of X remain bounded for all time. 
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To demonstrate the exponential divergences of X and ^ in the case of parametric 
resonance, recall from standard Floquet theory that there is a fundamental set of Jacobi 
fields that take the Bloch-Floquet forms Ji(T) = e^^ Ki{T) and J2(T) = e-^^i^slT), 
where Ki and K2 are periodic functions with the optical cycle being their common 
period [22]. Let = ciJi + C2J2, where ci and C2 are integration constants uniquely 
determined by the initial data for and vice verse. Thus diverges exponentially at 
the rate of the Floquet exponent /i for almost all initial data. Let < Ti < 1 be such 
that KiiTi) 7^ 0. Then for any integer > 0, 

J"(Ti + N)= ciKi(Ti)e'^(^i+^) + C2i^2(Ti)e-'^(^^+^). (5.60) 

Consequently, 

|cii^i(Ti)|e'^(^^+^) < I J^(Ti + N)\ + |c2i^2(ri)|e-'^(^^+^). (5.61) 

Hence if X{T) has exponential order a > for T > 0, then it follows from (15.571) and 
flCT]) that 

|cii^i(Ti)|e'^(^^+^) < M3e"(^i+^) + 102/^2(^1) le-'^^^^^^) 

or equivalently 

|cii^i(Ti)| < M3e("-'^)(^^+^) + |c2ir2(ri)|e-2^(^^+^) 

for some M3 > 0. If a < yU, this clearly leads to a contradiction as we pass to the 
limit N ^ 00 unless Ci = 0. Thus the oscillation center X must diverge exponentially 
faster than its associated Jacobi field . Indeed, the oscillatory component ^ must 
also diverge exponentially at a rate larger than than of . To see this, suppose ^ has 
exponential order /3 > for T > 0, then it follows from fl5.59p and fl5.6ip that 

|ciKi(Ti)|e^(^^+^) < M4e^(^^+^) + |c2K2(Ti)|e-'^(^^+^) 

or equivalently 

|ci/^i(Ti)| < M4e(^-'^)(^^+^) + |c2K2(Ti)|e-2'^(^^+^) 

for some M4 > 0. If /3 < /i, we also reach a contradiction by sending — )■ 00 unless 
ci = 0. 

To conclude this section, we note that the point of departure to validate the 
exponential divergences of dX/dT and d^/dT can be taken to be the equations 



^ = - V27rr/X(0) sin v^vrr/T + ^ 



cos Vo.TT'qT 

T=0 



d^ _ d^ 

It ~ It 



- 4:Ti^7f I COS V2Tir]{T - S) sin^ 2txS ■ i{S) dS, (5.62) 
Jo 

27rV / cos AttS ■ X{S)dS 
Jo 



T=0 

obtained by differentiating fl5.55p . Then precisely the same line of argument used to 
show the exponential divergences of X and ^ from 05.551) leads one from 05.621) to the 
conclusion that both dX/dT and d^/dT diverge exponentially faster than dJ^/dT for 
almost all initial data for in the case of a parametric resonance. 
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-1.084503 



Figure 1. The evolution of the Jacobi field and its associated Landau 
decomposition {X,(,) as well as the oscillation center X and its averaged and 
ponderomotive approximations X = Xa and Xp respectively for 10 optical cycles 
for different values of 77. The initial data are = 1 and dJ^ /dT = 0. When 77 ^ 1, 

« X « Xa ~ Xp for many cycles. When 77 = 1.2, the Jacobi field experiences 
a parametric resonance. In this case, both X and ^ diverge exponentially faster than 

. In consequence, neither Xa nor Xp serves as a reasonble approximation to X. 
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6. Conclusions 

In this paper, 

(i) we liave reformulated and analyzed the dynamical properties of the orbits of charged 
particles interacting with a generic E.M. field with planar symmetry, as governed by 
the Lorentz equations of motion, within a differential geometric framework whose 
geodesies suffer a mutual deviation in accordance with a curvature induced by the 
field intensity; 

(ii) we have demonstrated the integrability of the particle motion in a plane-wave field 
as a consequence of the vanishing of the curvature; 

(iii) we have indicated the methodology of examining the local stability of the orbit of 
a particle through the Jacobi field within such a geometrical formulation; and 

(iv) we have also showed the relevance of the geometrical formulation in discussing 
the domain of applicability of the ponderomotive oscillation center of a particle 
executing oscillatory motion in the E.M. field. 

In addition, by considering the motion of a charged particle in a linearly polarized 
standing wave field, we have shown that 

(a) different numerical subintervals of the impulse factor t] give rise to stable or unstable 
orbits as a consequence of the absence or the presence of parametric resonance of 
the Jacobi field; 

(b) in a sufficiently low impulsive E.M. field (77 ^ 1), a mathematically rigorous 
averaging principle can be applied to provide a precise meaning to the 
ponderomotive oscillation center of an orbit; 

(c) in the occurrence of parametric resonance, the applicability of the ponderomotive 
approximation breaks down. 

These are very interesting and promising results towards an in-depth differential 
geometric analysis for the dynamics of radiation-particle interaction. 
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